Molecular Dynamics Simulations of Hydroxyapatite Nanopores in Contact with Electrolyte Solutions: The Effect of Nanoconfinement and Solvated Ions on the Surface Reactivity and the Structural, Dynamical, and Vibrational Properties of Water

Hydroxyapatite, the main mineral phase of mammalian tooth enamel and bone, grows within nanoconfined environments and in contact with aqueous solutions that are rich in ions. Hydroxyapatite nanopores of different pore sizes (20 Å ≤ H ≤ 110 Å, where H is the size of the nanopore) in contact with liquid water and aqueous electrolyte solutions (CaCl2 (aq) and CaF2 (aq)) were investigated using molecular dynamics simulations to quantify the effect of nanoconfinement and solvated ions on the surface reactivity and the structural and dynamical properties of water. The combined effect of solution composition and nanoconfinement significantly slows the self-diffusion coefficient of water molecules compared with bulk liquid. Analysis of the pair and angular distribution functions, distribution of hydrogen bonds, velocity autocorrelation functions, and power spectra of water shows that solution composition and nanoconfinement in particular enhance the rigidity of the water hydrogen bonding network. Calculation of the water exchange events in the coordination of calcium ions reveals that the dynamics of water molecules at the HAP–solution interface decreases substantially with the degree of confinement. Ions in solution also reduce the water dynamics at the surface calcium sites. Together, these changes in the properties of water impart an overall rigidifying effect on the solvent network and reduce the reactivity at the hydroxyapatite-solution interface. Since the process of surface-cation-dehydration governs the kinetics of the reactions occurring at mineral surfaces, such as adsorption and crystal growth, this work shows how nanoconfinement and solvation environment influence the molecular-level events surrounding the crystallization of hydroxyapatite.